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A Quantum Field Theory Term Structure Model Applied to 

Hedging 

Belal E. Baaquie, Marakani Srikant, and Mitch Warachka 

Abstract 

A quantum field theory generalization, Baaquie M, of the Heath, Jarrow, and Morton 
(HJM) term structure model parsimoniously describes the evolution of imperfectly cor- 
related forward rates. Field theory also offers powerful computational tools to compute path 
integrals which naturally arise from all forward rate models. Specifically, incorporating field 
theory into the term structure facilitates hedge parameters that reduce to their finite fac- 
tor HJM counterparts under special correlation structures. Although investors are unable to 
perfectly hedge against an infinite number of term structure perturbations in a field theory 
model, empirical evidence using market data reveals the effectiveness of a low dimensional 
hedge portfolio. 



ca 1 Introduction 



■g : Applications of physics to hnance ate well known, and the application of qnantnm mechanics to the 



theory of option pricing is well known. Hence it is natural to utilize the formalism of quantum field 
theory to study the evolution of forward rates. Quantum field theory models of the term structure 
originated with Baaquie Q and the hedging properties of such models are analyzed in this paper. 
The intuition behind quantum field theory models of the term structure stems from allowing each 
forward rate maturity to both evolve randomly and be imperfectly correlated with every other 
maturity. As pointed out in Cohen and Jarrow |§ , this may also be accomplished by increasing the 
number of random factors in the original HJM towards infinity. However, the infinite number of 
factors in a field theory model are linked via a single function that governs the correlation between 
forward rate maturities. Thus, instead of estimating additional volatility functions in a multifactor 
HJM framework, one additional parameter is sufficient for a field theory model to instill imperfect 
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correlation between every forward rate maturity. As the correlation between forward rate maturities 
approaches unity, field theory models reduce to the standard one[| factor HJM model. Therefore, the 
fundamental difference between finite factor HJM and field theory models is the minimal structure 
the latter requires to instill imperfect correlation between forward rates. Section ^ demonstrates 
the challenge to reliably estimating volatility functions in a multifactor HJM model. However, it 
should be stressed that field theory models originate from the finite factor HJM framework with 
the Brownian motion(s) replaced by a field. As seen in the next section, crucial results such as the 
forward rate drift restriction are generalized but remain valid for an infinite factor process. 

The imperfect correlation between forward rates in a field theory model addresses the theoretical 
and practical challenges posed by their finite factor counterparts. Despite the enormous contribu- 
tion of the finite factor HJM methodology, an important dilemma is intrinsic to the finite factor 
framework. For a finite factor model to correctly match the movements of N forward rates, an 
iV factor model is required whose volatility parameters are difficult to estimate once N increases 
beyond one, Amin and Morton |T3| and section |j. However, field theory models only require an 
additional parameter to describe the correlation between forward rate maturities. Consequently, 
field theory models offer a parsimonious methodology to ensure the evolution of forward rates is 
consistent with the initial term structure without relying on a improbably large number of factors 
or eliminating the possibility of certain term structures, Bjork and Christensen More precisely, 
a field theory model does not permit a finite linear combination of forward rates to exactly describe 
the innovation in every other forward rate. From a practical perspective, an N factor term structure 
model implies the existence of an N dimensional basis for forward rates that allows N bonds with 
distinct maturities to be arbitrarily chosen when hedging. However, in the context of a field theory 
model, hedging performance depends on the maturity of the bonds in the hedge portfolio as illus- 
trated in subsection |4.1| . To summarize, the primary advantage of field theory models is their ability 
to generate a wider class of term structure innovations in a parsimonious manner. Consequently, 
hedging performance depends on the constituents of the hedge portfolio. This important issue is 
addressed empirically in section ^ after the derivation of hedge parameters in section |3|. As a final 
observation, stochastic volatility in a finite factor HJM model is implicitly a special case of field 
theory. 

Assigning independent random variables to each forward rate maturity has been studied by uti- 
lizing stochastic partial differential equations involving infinitely many variables. Past research that 
1 A field theory model may also converge to a multifactor HJM model as illustrated in Proposition [O. 
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instilled imperfect correlation between forward rates with generalized continuous random processes 
includes Kennedy [fT^I , Kennedy ]T3|, Goldstein ||, Santa Clara and Sornette [[RJ. The field theory 
approach developed by Baaquie is complimentary to past research since the expressions for all 
financial instruments are formally given as a functional integral. However, a significant advantage 
offered by field theory is the variety of computational algorithms available for applications involving 
pricing and hedging fixed income instruments. In fact, field theory was developed precisely to study 
problems involving infinitely many variables and the path integral described in the next section 
which serves as a generating function for forward curves is a natural tool to employ in applica- 
tions of term structure models. To date, there has been little empirical testing of previous field 
theory models which is understandable given their incomplete structure. The purpose of this paper 
is to demonstrate the practical advantages of utilizing field theory when hedging; specifically, the 
computational power offered by the path integral not found in previous research. 

Although field theory term structure models offer several improvements over finite factor mod- 
els, past research has not fully exploited all simplifying aspects of utilizing field theory. Specifically, 
implementation of field theory models when hedging bonds requires generalizing the concept of 
hedging. Since a field theory model may be viewed as an infinite factor model, perfectly hedging 
a bond portfolio is impossible?] although computing hedge parameters remains feasible. Further- 
more, empirical evidence in section f| finds the correlation between forward rate maturities exerts a 
significant impact on the hedging of bonds and allows a low dimensional basis to effectively hedge 
term structure risk. It is important to emphasize that field theory facilitates rather than inhibits 
the implementation of term structure models as seen when the path integral for forward rates is 
introduced. The path integral is a powerful theoretical tool for computing expectations involving 
forward rates required for pricing contingent claims such as futures contracts in a straightforward 
manner. As expected, closed form solutions for fixed income contingent claims and hedge parame- 



ters reduce to standard "textbook" solutions such as those found in Jarrow and Turnbull |l 1] when 
the correlation between forward rates approaches unity. In summary, this paper elaborates on an 
implementable field theory term structure model that addresses the limitations inherent in finite 
factor term structure models. 

This paper is organized as follows. Section two briefly summarizes the field theory model with 
a detailed presentation found in the appendix. The third and fourth sections cover the theoretical 
and empirical aspects of hedging in a field theory model. The conclusion is left to section five. 



2 A measure valued trading strategy as in Bjork, Kabanov, and Runggaldier || provides one alternative. 
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2 Field Theory 



The field theory model underlying the hedging results in section three was developed by Baaquie 
and calibrated empirically in Baaquie and Marakani J3|. As alluded to earlier, the fundamental 
difference between the model presented in this paper and the original finite factor HJM model stems 
from the use of a infinite dimensional random process whose second argument admits correlation 
between maturities. 

A Lagrangian is introduced to describe the field. The Lagrangian has the advantage over 
Brownian motion of being able to control fluctuations in the field, hence forward rates, with respect 



to maturity through the addition of a maturity dependent gradient as detailed in Definition |27[ 



The action of the field integrates the Lagrangian over time and when exponentiated and normal- 
ized serves as the probability distribution for forward rate curves. The propagator measures the 
correlation in the field and captures the effect the field at time t and maturity x has on maturity x' 
at time t' . In the one factor HJM model, the propagator equals one which allows the quick recovery 
of one factor HJM results. Previous research by Kennedy []I2"| |, Kennedy and Goldstein Hj has 
begun with the propagator or "correlation" function for the field instead of deriving this quantity 
from the Lagrangian. More importantly, the Lagrangian and its associated action generate a path 
integral that facilitates the solution of contingent claims and hedge parameters. However, previous 
term structure models have not defined the Lagrangian and are therefore unable to utilize the path 
integral in their applications. The Feynman path integral, path integral in short, is a fundamen- 
tal quantity that provides a generating function for forward rate curves and is formally introduced 
in Definition |A.1| of the appendix. Although crucial for pricing and hedging, the path integral has 
not appeared in previous term structure models with generalized continuous random processes. 

Remark 2.1 Notation 

Let to denote the current time and T the set of forward rate maturities with to < T. The upper bound 
on the forward rate maturities is the constant Tfr which constrains the forward rate maturities T 
to lie within the interval [t ,t + T FR ]. 

To illustrate the field theory approach, the original finite factor HJM model is derived using field 
theory principles in appendix [A]. In the case of a one factor model, the derivation does not involve 
the propagator as the propagator is identically one when forward rates are perfectly correlated. 
However, the propagator is non trivial for field theory models as it governs the imperfect correlation 
between forward rate maturities. Let A(t, x) be a two dimensional field driving the evolution of 
forward rates f(t,x) through time. Following Baaquie the Lagrangian of the field is defined as 



Definition 2.1 Lagrangian 
The Lagrangian of the field equals 




2T FR I ' fx 2 \ dx 
Definition [2.1| is not unique, other Lagrangians exist and would imply different propagators. 



However, the Lagrangian in Definition 2.1 is sufficient to explain the contribution of field theory 



while quickly reducing to the one factor HJM framework. Observe the presence of a gradient with 

dACt x) 

respect to maturity — — that controls field fluctuations in the direction of the forward rate 

ox 

maturity. The constant^ /i term measures the strength of the fluctuations in the maturity direction. 
The Lagrangian in Definition |2.1| implies the field is continuous, Gaussian, and Markovian. Forward 
rates involving the field are expressed below where the drift and volatility functions satisfy the usual 
regularity conditions. 



= a{t,x) + a{t,x)A(t,x) (2) 

The forward rate process in equation incorporates existing term structure research on Brown- 
ian sheets, stochastic strings, etc that have been used in previous continuous term structure models. 
Note that equation (§) is easily generalized to the K factor case by introducing K independent and 
identical fields Ai(t, x). Forward rates could then be defined as 

df_t± x) = ^ x) + ^2 x)Ai(t, x) (3) 
i=i 

However, a multifactor HJM model can be reproduced without introducing multiple fields. In 
fact, under specific correlation functions, the field theory model reduces to a multifactor HJM model 
without any additional fields to proxy for additional Brownian motions. 

Proposition 2.1 Lagrangian of Multifactor HJM 

The Lagrangian describing the random process of a K-f actor HJM model is given by 



C[A] = --A(t,x)G- 1 (t,x,x')A(t,x') 



'Other functional forms are possible but a constant is chosen for simplicity. 
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where 



df(t.x) , . , . 

K m ' = a(t,x)+A(t,x) 

and x, x') denotes the inverse of the function 

K 

G(t,x,x') = o-j(t, x)o-j(t, x) 

i=l 

The above proposition is an interesting academic exercises to illustrate the parallel between 
field theory and traditional multifactor HJM models. However, multifactor HJM models have 
the disadvantages described earlier in the introduction associated with a finite dimensional basis. 
Therefore, this approach is not pursued in later empirical work. In addition, it is possible for 
forward rates to be perfectly correlated within a segment of the forward rate curve but imperfectly 
correlated with forward rates in other segments. For example, one could designate short, medium, 
and long maturities of the forward rate curve. This situation is not identical to the multifactor 
HJM model but justifies certain market practices that distinguish between short, medium, and 
long term durations when hedging. However, more complicated correlation functions would be 
required; compromising model parsimony and reintroducing the same conceptual problems of finite 
factor models. Furthermore, there is little economic intuition to justify why the correlation between 
forward rates should be discontinuous. Therefore, this approach is also not considered in later 
empirical work. 

2.1 Propagator 

The propagator is an important quantity that accounts for the correlation between forward rates 
in a parsimonious manner. The propagator D(x,x';t,TFR) corresponding to the Lagrangian in 



Definition [2.1| is given by the following lemma where #(•) denotes a Heavyside function. 



Lemma 2.1 Evaluation of Propagator 
The propagator equals 
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D(x, x'\ t, T, 



FR 



FR) 



sinh fJ,(T FR — t) sinh(/ir')6 l (r — r') 



siah(fj,T FR ) 
+ sinh fi(T FR — r') sinh(/xr)^(r' — r) 



+ o , rr -{2 COSh /i \T — I COSh U ( T — ] 

2cosh 2 (^)l *\ 2 J 2 y 

+ sinh(/ir) sinh( / ur / ) + sinh fi(T FR — t) sinh ft(T FR — r') | 

where t = x — t and r 1 — x' — t both represent time to maturities. 



Lemma [2.1| is proved by evaluating the expectation E[A(t,x),A(t',x')\. The computations are 
tedious and contained in Baaquie JI| but are well known in physics and described in common 
references such as Zinn- Justin |p7| . The propagator of Goldstein [0 is seen as a special case of 



Lemma 2A defined on the infinite domain —00 < x, x' < 00 rather than the finite domain t < 
x,x' < t + TpR. Hence, the propagator in Lemma |2.1| converges to the propagator of Goldstein 
as the time domain expands from a compact set to the real line. The effort in solving for the 
propagator on the finite domain is justified as it allows covariances near the spot rate f(t, t) to differ 
from those over longer maturities. Hence, a potentially important boundary condition defined by 
the spot rate is not ignored. 

Observe that the propagator D(x, x'; t, Tp^) in Lemma [2.1| only depends on the variables r and 
t' as well as the correlation parameter /i which implies that the propagator is time invariant. This 
important property facilitates empirical estimation in section [| when the propagator is calibrated 
to market data. To understand the significance of the propagator, note that the correlator of the 
field A(t, x) for to < t, t' < to + Tfr is given by 



E[A(t,x)A(t',x')} = 5(t-t')D(x,x';t,T FR ) (4) 

In other words, the propagator measures the effect the value of the field A(t, x) has on A(t' , x'); 
its value at another maturity x' at another point in time. Although D(x, x'\ t, T FR ) is complicated in 
appearance, it collapses to one when \x equals zero as fluctuations in the x direction are constrained 
to be perfectly correlated. It is important to emphasize that fi does not measure the correlation 
between forward rates. Instead, the propagator solved for in terms of fi fulfills this role. 

Remark 2.2 Propagator, Covariances, and Correlations 

The propagator D(x, x'; t, T FR ) serves as the covariance function for the field while 

a(t, x)D(x, x'] t, T FR )a(t, x ) 
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serves as the covariance function for forward rates innovations. Hence, the above quantity is repeat- 
edly found in hedging and pricing formulae presented in the next section. The correlation functions 
for the field and forward rate innovations are identical as the volatility functions a(t, ■) are elim- 
inated after normalization. This common correlation function is critical for estimating the field 
theory model and equals equation of section [|. 

The following tablef] summarizes the important terms in both the original HJM model and its 
extended field theory version. 



Quantity 


Finite Factor HJM 


Field Theory 


Lagrangian 


-\W\t) 


i J 

2T FR 




Propagator 


1 


D(x,x';t,T FR ) 


Path Integral 


exp{i/~dtf»(t)} 


exp 1 1 dt jl +TpR dxdx'J(t, x)D(x, x'\ t, T FR )J(t, x') j 



As expected, the HJM drift restriction is generalized in the context of a field theory term struc- 
ture model. However, producing the drift restriction follows from the original HJM methodology 
as the discounted bond price evolves as a martingale under the risk neutral measure to ensure no 
arbitrage. Under the risk neutral measure, the bond price is written as 



P(t ,T) = E\ 



[to,U 



•J£r®dt 



P(U,T) 



VAe~ dtm) e~ £ dxf{t *' x) 



(5) 



where j T>A represents an integral over all possible field paths in the domain J** dt dx. The 
notation -E[t ,tj[S'] denotes the expected value under the risk neutral measure of the stochastic 
variable S over the time interval [to>£*]- Equation (||) serves as the foundation for computing the 
forward rate drift restriction stated in the next proposition and proved in appendix |A|. 

Proposition 2.2 Drift Restriction 

The field theory generalization of the HJM drift restriction equals 



a(t, x) 



a(t,x) J dx'D(x,x';t,T FR )a(t,x') 



4 The j(t) and J(t, x) functions found in the definition of the path integral are "source" functions used to compute 
the moments of forward curves (see appendix |Aj) . They do not appear in the solution of contingent claims or hedge 
parameters. 
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As expected, with \i equal to zero the result of Proposition ^]2] reduces to 

a(t, x) — a(t, x) I dx'a(t, x') 



and the one factor HJM drift restriction is recovered. The next section considers the problem of 
hedging in the context of a field theory model using either bonds or futures contracts on bonds. 

3 Pricing and Hedging in Field Theory Models 

Hedging a zero coupon bond denoted P(t, T) using other zero coupon bonds is accomplished by 
minimizing the residual variance of the hedged portfolio. The hedged portfolio II(t) is represented 
as 

N 

i=i 

where Aj denotes the amount of the i th bond P(t, Tj) included in the hedged portfolio. Note 
the bonds P(t,T) and P(t, 2$) are determined by observing their market values at time t. It is 
the instantaneous change in the portfolio value that is stochastic. Therefore, the volatility of this 
change is computed to ascertain the efficacy of the hedge portfolio. 

For starters, consider the variance of an individual bond in the field theory model. The definition 
P(t, T) = exp (— J T dxf(t, x)) for zero coupon bond prices implies that 



dP(t,T) 
P(t,T) 



f(t,t)dt 



dxdf(t, x) 



r 



r(t) — / dxa(t,x) — / dxa(t, x)A(t, x) dt 



and E 



dP(t,T) 

p(t,t) 



r(t) — dxa(t,x)) dt since E[A(t,x)\ = 0. Therefore 



dP(t,T) 
P(t,T) 



E 



dP(t,T) 
[P(t,T) 



—dt / dxa(t, x) A(t, x) 



(6) 



Squaring this expression and invoking the result that E[A(t,x)A(t,x r )] = S(0)D(x,x r ;t,TFR) 
D< " x ' x f t ' TFR ' > results in the instantaneous bond price variance 



Var[dP(t,T)} = dtP 2 (t,T) f dx [ dx'a(t,x)D(x,x';t,T FR )a(t,x') 

Jt Jt 



(7) 
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As an intermediate step, the instantaneous variance of a bond portfolio is considered. For a 
portfolio of bonds, fl(t) = ^ i=1 AjP(i, Tj), the following results follow directly 



dfi{t) - E[dU{t)} = -dt V AiP{t, Ti) j dxa(t, x)A(t, x ) 

„-_i Jt 



i=l 

and 

N N 



r 1 i r 1 ] 

Var[dfl(t)] = dt^2^2A i A j P(t,T l )P(t,T j ) dx dxa(t,x)D(x,x';t,T FR )a(t,x') (9) 
i=l j=i Jt Jt 

The (residual) variance of the hedged portfolio 

N 

i=i 

may now be computed in a straightforward manner. For notational simplicity, the bonds P(t, TJ) 
(being used to hedge the original bond) and P(t, T) are denoted Pi and P respectively. Equation 
(H) implies the hedged portfolio's variance equals the final result shown below 

P 2 dx dx'a(t, x)a(t, x')D(x, x';t,T FR ) 

N rT rTi 

+2P s ^A i P i J dx J dx'a(t, x)a(t, x')D(x, x';t,TFn) 

N N Ti T . 

+ / ] / ] AjAjPjPj / dx dx'a(t, x)a(t, x')D(x,x';t,TFn) 

i=l j=l J* J* 

Observe that the residual variance depends on the correlation between forward rates described 
by the propagator. Ultimately, the effectiveness of the hedge portfolio is an empirical question 
since perfect hedging is not possible without shorting the original bond. This empirical question is 
addressed in section |] when the propagator is calibrated to market data. Minimizing the residual 
variance in equation ([II]) with respect to the hedge parameters A$ is an application of standard 
calculus. The following notation is introduced for simplicity. 

Definition 3.1 

rT rT 

Li = PPi I dx I dx'a(t,x)a(t,x')D(x,x';t,T FR ) 
Jt Jt 

Mij = PiPj / dx dx a(t, x)a(t,x')D(x, x ;t,T FR ) 
Jt Jt 
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Definition 



3.1 



allows the residual variance in equation ([[(]) to be succinctly expressed as 



N N N 



P 2 j dx j dx'a(t,x)a{t,x')D(x,x']t,T FR ) + 2^A i L i + ^2^TA i A j M ij (11) 

Jt J t i=l i=l j=l 

Hedge parameters A, that minimize the residual variance in equation ( p~T|) are the focus of the 
next theorem. 

Theorem 3.1 Hedge Parameter for Bond 
Hedge parameters in the field theory model equal 



N 

and represent the optimal amounts of P(t, Tj) to include in the hedge portfolio when hedging 
P{t,T). 

Theorem [Tl] is proved by differentiating equation (|ll|) with respect to Aj and subsequently 
solving for Aj. Corollary |3.1| below is proved by substituting the result of Theorem |3.1| into equation 

fITTD. 



Corollary 3.1 Residual Variance 

The variance of the hedged portfolio equals 



V = P 2 j dx j dx'o-(t,x)a(t,x')D(x,x';t,T FR ) -^^LiMr^Lj 

Jt Jt i=l j=l 

which declines monotonically as N increases. 



The residual variance in Corollary |3.1| enables the effectiveness of the hedge portfolio to be 
evaluated. Therefore, Corollary |3.1| is the basis for studying the impact of including different bonds 
in the hedge portfolio as illustrated in subsection [O . For N — 1, the hedge parameter in Theorem 



3.1 reduces to 



A _ _P_ ( dx J? 1 dx'ajt, x)a(t, x')D(x, x'; t, T FR ) ^ 
Pi \ J^ 1 dx J^ 1 dx'a(t, x)a(t, x')D(x, x'\ t, T FR ) 
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To obtain the HJM limit, constrain the propagator to equal one. The hedge parameter in 
equation ( |I2] ) then reduces to 



Ai 



P 



^ dx J*. Tl dx'a(t, x)a(t, x 



V 



jj 1 dxa(t. 



x 



I 



P I dxa(t,x) 
Pi \ J^ 1 dxa(t, x) 



(13) 



The popular exponential volatility function a(t, T) = ae~ x( - T ~^ allows a comparison between 
our field theory solutions and previous research. Under the assumption of exponential volatility, 
equation (|13|) becomes 



Ax 



P / i _ e -A(T-t) 
T x VI -e-^-') 



(14) 



Equation (14) coincides with the ratio of hedge parameters found as equation 16.13 of Jarrow 



and Turnbull [FTTFI . In terms of their notation 



Ax 



P(t,T) fX(t,T) 



P(t,T0 V*M) 
For emphasis, the following equation holds in a one factor HJM model 

d[P(t,T) + A 1 P(t,T 1 )} 



(15) 



dr{t) 







which is verified using equation ([15]) and results found on pages 494-495 of Jarrow and Turnbull 

m 



d[P(t 1 T) + A 1 P(t,T 1 ) 
dr(t) 



-P(t, T)X(t, T) - A 1 P(t, Tt)X(t, Ti) 
-P(t, T)X(t, T) + P(t, T)X(t, T) = 



When Ti = T, the hedge parameter equals minus one. Economically, this fact states that the 
best strategy to hedge a bond is to short a bond of the same maturity. This trivial approach 
reduces the residual variance in equation (O) to zero as Ai = — 1 and P = Pi implies L\ = M n . 



Empirical results for nontrivial hedging strategies are found in subsection |44] after the propagator 
is calibrated. 
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3.1 Futures Pricing 

As this paper is primarily concerned with hedging a bond portfolio, futures prices are derived as 
they are commonly used for hedging bonds given their liquidity. 

Proposition 3.1 Futures Price 

The futures price J 7 {to, £*, T) is given by 



F{t ,U,T) = E [tQ>u] [P{U,T)\ 



= F(t ,U,T)exp{n T (t ,U,T)} (16) 
where F(to,t*,T) represents the forward price for the same contract, F(t ,t*,T) = and 



rU ft, rT 

^V(*o> t*i T) = — / dt dxa(t,x) / dx D(x,x ;t,T FR )a(t, x') (17) 
Jtn Jt Ju 



r-t, rt, rT 

dt / dxcr(t,x) 

1 1 Jt Ju 

Details of the proof are found in appendix [A]. Observe that for fi = 0, equation ( |TT| ) collapses to 

pt, i>t, pT 

VLjr(tQ,t*,T) = — / dt I dxcr(t,x) / dx'a(t,x') (18) 
Jt Jt Jt, 

which is equivalent to the one factor HJM model. For the one factor HJM model with exponential 
volatility, equation fllTf) becomes 



^(£0, t*, T) 



a 



2A 3 



(1 _ e -A(T-t»)) ^ _ e 



-A(t„,-t )) 2 



which coincides with equation 16.23 of Jarrow and Turnbull [TT[. Observe that the propagator 
modifies the product of the volatility functions with serving as an additional model parameter. 
Prices for call options, put options, caps, and floors proceed along similar lines with an identical 
modification of the volatility functions but their solutions are omitted for brevity. Formulae for 
these contingent claims are given in Baaquie M where their reduction to the closed form solutions 



of Jarrow and Turnbull |ll]] is also presented. 
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3.2 Hedging Bonds with Futures Contracts 

The material in the previous subsection allows the hedging properties of futures contracts on bonds 
to be studied. Proceeding as before, the appropriate hedge parameters for futures contracts expiring 



in one year are computed. Proposition |3J] expresses the futures price T(t, t*, T) in terms of the 
forward price p^'f^ = e~ ^* dx ^ t,x ^ and the deterministic quantity £V(£, t*,T) found in equation 
fllTD- The dynamics of the futures price dT(t, £*, T) is given by 



^UT) = dQAt > U ' T) ~ f dxdf{t ' X) (19) 



which implies 



mt,t„T) - E[ dH t,t^n = _ dt f T dMtx)A(ttX) (20) 



u 



F(t,t*,T) 

Squaring both sides leads to the instantaneous variance of the futures price 



Var[cLF(t,t*,T)] = dtF 2 (t,t*,T) J dx J dx'cr(t,x)D(x,x)a(t,x') (21) 



fT i-T 

dx I 

The following definition updates Definition p.l| in the context of futures contracts. 



Definition 3.2 Futures Contracts 

Let Ti denote the futures price !F(t,t*,Ti) of a contract expiring at time t* on a zero coupon bond 
maturing at time Ti . The hedged portfolio in terms of the futures contract is given by 

N 

n(t) = p+j2 A iFi 

where Ti represent observed market prices. For notational simplicity, define the following terms 



f f 

Li = PTi I dx I dx , a(t,x)D(x,x';t,TFji)a(t,x) 

dx / dx'o-(t,x)D(x,x']t*iTFR)°~{t-, x ') 
The hedge parameters and the residual variance when futures contracts are used as the under- 



lying hedging instruments have identical expressions to those in Theorem 3.1 and Corollary |3.1| but 



are based on Definition 372. Computations parallel those in the beginning of this section. 
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Corollary 3.2 Hedge Parameters and Residual Variance using Futures 

Hedge parameters for a futures contract that expires at time on a zero coupon bond that matures 
at time Ti equals 



N 



while the variance of the hedged portfolio equals 



r-T pT n n 

V = P 2 dx dx'a(t,x)a(t,x')D(x,x';t,T FR ) -^^LiMr^Lj 



i=i j=i 



for Li and Mij in Definition 3.2 



Proof follows directly from previous work. 

4 Empirical Estimation of Field Theory Models 

This section illustrates the significance of correlation between forward rate maturities and numeri- 
cally estimates its impact on hedging performance. The volatility function o~(t, x) and the /x param- 
eter were previously calibrated non parametrically from market data in Baaquie and Srikant ||. 
The data used in the following empirical tests was generously provided by Jean-Philippe Bouchaud 
of Science and Finance. The data consists of daily closing prices for quarterly Eurodollar futures 
contracts with a maximum maturity of 7.25 years as described in Bouchaud, Sagna, Cont, and 
El-Karoui [HI as well as Bouchaud and Matacz HI . 



As previous empirical research into estimating HJM models has found, Amin and Morton [15|, 
multifactor HJM models are difficult to estimate. This concept is revealed by considering the 
difference between a one and two factor HJM model. Let r = x — t and r' — x' — t once again 
represent time to maturity and assume the correlation function is time invariant by depending solely 
on r and r'. The correlation between forward rate maturities C(x,x') may be expressed in terms 
of the normalized propagator R as 



C(x,x') 



D(x, x'; t) 
D(x, X] t)D(x', x'; t) 



5 The value of Tfr was set to oo, rescaling by ^i— preserved the Lagrangian in Definition 
propagator. 



2.1 



(22) 

and its associated 
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To clarify, the propagator (as emphasized after Lemma |2.1|) and hence the correlation in equa- 
tion fl2^D is time homogeneous. Once the propagator is determined, the correlation between forward 
rates follows immediately and vice versa. For example, when the propagator equals one, the corre- 
lation between forward rates is also identically one. However, for the two factor HJM model, the 
normalized propagator is given by 



l+g(r)g(r') 



(23) 



Consequently, the correlation function above in equation (El) depends on the 



for 9(r) = 

ratio of volatility functions. Therefore, for a given correlation structure, obtaining reliable volatility 
function estimates is a challenge as it is difficult to disentangle one volatility function from another. 
This practical disadvantage is overcome by field theory models. 

Concerning the estimation of the field theory model, the market volatility function is used 
throughout the remainder of section f|. This volatility function was estimated directly from market 
data as the variance of term structure innovations for each maturity and is graphed below. 



Sigma estimated from data 



0.06 



0.055 



0.05 



0.045 



0.04 



0.035 




3 4 5 

Time to maturity/year 

Figure 1: Implied Volatility Function using Market Data 

The corresponding implied correlation parameter /j, was estimated as .06 (annualized). Esti- 
mation of \i was accomplished after estimating the correlation between forward rate innovations 
over the 7.25 year horizon. The fi parameter was then found by minimizing the root mean square 
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difference between the theoretical correlation function in equation (p2) which contains \x via the 
definition of the propagator in Lemma 271 and the empirical correlation function generated from the 
data aggregated over the sample period. The minimization was carried out using the Levenberg- 
Marquardt method. Stability analysis performed by considering different subsections of the data 
indicated that this estimate was robust with an error of at most 0.01. The propagator itself is 
graphed below. 



Correlation between different maturities 



line 1 



correlation 



1 

0.95 
0.9 

0.85 
0.8 

0.75 




maturity/year 



maturity/year 



Figure 2: Propagator Implied by \i 

It should be emphasized that the propagator may also be estimated non parametrically from the 
correlation found in market data without any specified functional form, as the volatility function 
was estimated. This approach preserves the closed form solutions for hedge parameters and futures 
contracts illustrated in the previous section. However, the original finite factor HJM model cannot 
accommodate an empirically determined propagator since it is automatically fixed once the HJM 
volatility functions are specified. 

4.1 Hedging Error 

The reduction in variance achievable by hedging a five year zero coupon bond with other zero coupon 
bonds and futures contracts is the focus of this subsection. The residual variances for one and two 
bond hedge portfolios are shown in figures |3| and [|. The parabolic nature of the residual variance is 
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because n is constant. A more complicated function would produce residual variances that do not 
deviate monotonically as the maturity of the underlying and the hedge portfolio increases although 
the graphs appeal to our economic intuition which suggests that correlation between forward rates 
decreases monotonically as the distance between them increases as shown in figure [| Observe that 
the residual variance drops to zero when the same bond is used to hedge itself; eliminating the 
original position in the process. The corresponding hedge ratios are shown in figure |^. 
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Residual variance when hedging 5 year zero coupon bonds 
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Figure 3: Residual variance for five year bond versus bond maturity used to hedge. 



Residual variance when hedging 5 year zero coupon bonds 



line 1 




3 Time to maturity/year 



Time to maturity/year 



Figure 4: Residual variance for five year bond versus two bond maturities used to hedge. 
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Hedge ratio when hedging a five year zero coupon bonds 
I 1 1 1 1 1 1 — 



-5 - 



-10 - 



-15 - 



-20 - 




_25 I i i i i i i i I 

012345678 

Time to maturity/year 

Figure 5: Hedge ratios for five year bond. 

Numerical tests to determine the efficiency of hedging using futures contracts that expire in one 
year are conducted by calculating the residual variance when hedging a five year zero coupon bond 
with a futures contract expiring in one year on various zero coupon bonds. The residual variance 
is shown in figure ^. Observe that the zero coupon bond is best hedged by selling futures contracts 
on 4.5 year bonds which is explained by the fact that the futures contract only depends on the 
variation in forward rate curve from £* to T but the zero coupon bond depends on the variation 
of the forward rate curve from t to T. Hence, a shorter underlying bond maturity is chosen for 
the futures contract to compensate for the forward rate curve from t to t*. A similar result is seen 
when hedging with two futures contracts both expiring in one year. In this case, optimal hedging 
is obtained when futures contracts on the same bond as well as one on a bond with the minimum 
possible maturity (1.25 years) are shorted. The use of a futures contract on a short maturity bond 
is consistent with the high volatility of short maturity forward rates as displayed in figure [I]. The 
optimal futures contracts to include in the hedge portfolio are shown in table |I] when hedging a 5 
year bond. 
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Residual variance when hedging with one futures contract 
0.0014 | 1 1 1 1 1 



0.0012 - 



0.001 - 



0.0008 - 



0.0006 - 



0.0004 - 



0.0002 - 





1 2 3 4 5 6 7 8 

Time to maturity/year 

Figure 6: Residual variance for five year bond hedged with a one year futures contract on a T 
maturity bond. 



Number 


Futures Contracts (Hedge Ratio) 


Residual Variance 





none 


0.03046 


1 


4.5 years (-1.288) 


8.81 x 10" 5 


2 


5 years (-0.9347), 1.25 years (-2.72497) 


2.63 x 10" 5 



Table 1: Residual variance and hedge ratios for a five year bond hedged with one year futures 
contracts. 

The variances recorded in Table |l| correspond to daily dollar denominated fluctuations. 
4.2 Current and Future Research 

The field theory model has been used to value call and put options as well as caps and floors in 
Baaquie ffl. Hedge parameters for these instruments is currently under investigation. Obtaining 
fixed income options data to investigate the significance of correlated forward rates on option prices 
as well as its impact on their hedge parameters would also be valuable. An enhanced model with 
stochastic volatility where innovations in the term structure result from the product of two fields 
has been developed by Baaquie ||. 
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5 Conclusion 



Field theory models address the theoretical limitations of finite factor term structure models by 
allowing imperfect correlation between every forward rate maturity. The field theory model offers 
computationally expedient hedge parameters for fixed income derivatives and provides a methodol- 
ogy to answer crucial questions concerning the number and maturity of bonds to include in a hedge 
portfolio. Furthermore, field theory models are able to incorporate correlation between forward rate 
maturities in a parsimonious manner that is well suited to empirical implementation. Empirical 
evidence revealed the significant impact that correlation between forward rate maturities has on 
hedging performance. Despite the infinite dimensional nature of the field, it is shown that a low di- 
mensional hedge portfolio effectively hedges interest rate risk by exploiting the correlation between 
forward rates. Therefore, field theory models address the theoretical dilemmas of finite factor term 
structure models and offer a practical alternative to finite factor models. 



A Details of Field Theory Model 

This appendix briefly reviews the results contained in Baaquie [[TJ. The formalism of quantum 
mechanics is based on conventional mathematics of partial differential equations and functional 



analysis as detailed in Zinn- Justin HPT |. The mathematical tools underlying quantum field theory 
have no counterpart in traditional stochastic calculus although Gaussian fields are equivalent to an 
infinite collection of stochastic processes. The Lagrangian has the advantage over Brownian motion 
of being able to control fluctuations in the field, hence forward rates, with respect to maturity 
through the addition of a maturity dependent gradient as detailed in equation ([H]). The action 



integrates the Lagrangian over time and when exponentiated and appropriately normalized yields a 
probability distribution function for the forward rate curves. The propagator measures the correla- 
tion in the field and captures the effect the field at time t and maturity x has on the maturity x' at 
time t' . The Feynman path integral serves as the generating function for forward rate curves. The 
path integral is obtained by integrating the exponential of the action over all possible evolutions of 
the forward rate curve. 



A.l Restatement of HJM 

Before presenting the field theory of the forward rates, we briefly restate the HJM model in terms 
of its original formulation but using the notation of path integrals. For simplicity, consider a one 
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factor HJM model of the forward rate curve whose evolution is generated by 



df ^ X) = a(t, x) + a(t, x)W(t) (24) 

where a(t, x) and a(t, x) are the drift and volatility of forward rates. For every value of time t, 
the stochastic variable W(t) is an independent Gaussian random variable with the property that 

E[W{t)W{t')) = 8{t - f) (25) 

in contrast to equation (§) involving the field. More conventionally, the HJM model is written 
as df(t,x) = a(t,x)dt + a(t, x)dZ(t) where Z(t) represents a Brownian motion process. Hence, the 
Gaussian process W(t) equals ^jr-- To derive the covariance function implied for W(t), the corre- 
lation function for Z(t) must be twice differentiated. Therefore, {W(t)W(t')) = -^j(Z(t)Z(t')) = 

d 2 min(t/) _ dB{t-t') 



;jl - 5(t — t') as seen in equation (p5|). The Lagrangian of the Gaussian process is 
defined as 



Cm = -\w\t) (26) 

The Gaussian process may be illustrated by discretizing time into a discrete lattice of spacing 
e, and t — > me, with m = 1,2....M, with W(t) — > W(m). The probability measure underlying 
Gaussian paths for t\ < t < t 2 is given by 



M M 



V[W] = J\^ W2{m) = Y[e eC[w] = e^^c[w] (27 ) 

m=l Tn=l 

M j— r+oo 

m=l v " /_ °° 



/ 



The term PfW] represents the probability of a path traced out by W(t). For purposes of rigor, 
the continuum notation represents taking the continuum limit of the discrete multiple integrals 

/+oo 
dx(t) has a rigorous, 

measure theoretic, definition as the integration over all continuous, but nowhere differentiable, paths 
running between points x(t\) and xfa)- For t\ < t < the probability distribution function for 
the paths of the Gaussian process equals 
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V[WMM=™-V{S[WMM} (28) 

where the action S[W, ti, t 2 ] f° r the Gaussian process is given by the limit of the exponent in 
equation ( p7f) as e — > 

S , [W,ti,t 2 ] = ~jT* a w 2 (t)dt (29) 

The term J PVV denotes path integration over all the random variables W(t) which appear in 
the problem. A path integral approach to the HJM model has been discussed in Chiarella and 
El-Hassan [0 although the action derived is different than the one given above since a different 
set of variables were involved. A formula for the generating function of forward rates driven by a 
Gaussian process is given by the path integral 



Z[jMM = J VWe^ dtmw{t) e s ^^ 



ii 



\\ I fin<lt} (30) 



This path integral is crucial for applications involving the pricing of derivatives as demonstrated 
in section |3|. 

A. 2 Field Theory Model 

The Lagrangian of the field equals 



Since the Lagrangian is quadratic, the resulting field is Gaussian and the existence of a Hamil- 
tonian for the Lagrangian implies the field is Markovian. Forward rates are expressed as 



= a(t,x) + a{t,x)A{t,x) (32) 
The action of S[A] integrates the Lagrangian over time to yield 
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S[A] = / dt dxC[A] 

J to Jt 

rt+T FR ( 

dt dxl A 2 (t, x) + ^[ ) K (33) 



l r°° , f t+TpR ,1,2, , l /<M(t,x^ 2 



Consider the one factor case 



= / A{t,x)dx (34) 



by assuming the correlation fi is zero. This process reduces the action found in equation fl33|) to 
the following action 



I poo ft+T FR 

S[A] = -— / W\t)dt / dx 



Ffi J io Jt 

oo 

VT 2 (f)df (35) 

*0 



which is identical to equation (|29|). If one thinks of the field A(to,x) at some instant to as the 
position of a string, then the one factor HJM model constrains the string to be rigid. The action 
S[A] given in (|3^) allows every maturity x in A(t ,x) to fluctuate as a string with string rigidity 
equal to -\. For [i — 0, the string is infinitely rigid or perfectly correlated. A normalizing constant 
Z is the result of integrating the paths of the field against the probability of each path. 



VAe s[A] (36) 

where the notation J T>A represents an integral over all possible field paths contained in the 
domain. The moment generating functional for the field is given by the Feynman path integral, 



Zinn- Justin HIT 



Definition A.l Path Integral 

The path integral over all possible paths of the field, weighted by their probability, equals 

Z [J] = iy VAe! JA e s ^ 
where Z is defined in equation ($0). 
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In the context of the present term structure model, the path integral equals 



i p C r<x> pt+T FR -\ 

Z[J] = -JVAexp^J^ dt dxJ(t,x)A(t,x)je slA] (37) 

The above formula is interpreted as follows, the Z term operates as a normalizing constant while 

S[A] 

^-g— represents the probability distribution function for each random path generated by the field 
A(t,x). The integrand 

f roo rt+T FR >. 

exp < / dt I dxJ(t,x)A(t,x) > 

contains an external source^] function J(t, x) coupled to the field A(t, x) that is operated on to 
produce moments of the distribution. The path integral Z\J\ is evaluated explicitly in Baaquie 
and simplifies equation (|37|) to 



f 1 f°° pt+TFR 

Z[J] = exp|- J dt J dxdx'J(t,x)D(x,x';t,TFR)J(t,x'] 



(3? 



The path integral is the basis for all applications of term structure modeling such as the pricing 
and hedging of fixed income contingent claims since it represents the generating function for forward 
rate curves. 

A. 3 Proof of Proposition |ZTZ] 

Starting from equation (H) 

P(t ,T) = JvAe-fiWMe-Z* 1 '^ 
and using the identity P(to,T) = e"^o dxj ^ to ' x ^ implies 

1 = J VAe'^o 
1 = e" 

and isolates the drift in the exponent as 



- /** dt J t T dxa(t,x)+a(t,x)A(t,x) 



- ff* dt J t T dxa(t,x)+\ /** dt dx dx' <j(t,x)D(x,x' ;t,T FR )a(t,x') 



6 To make an analogy with univariate generating functions, the function J(t, x) serves a similar role to t in the 
moment generating function e t J - t +k< 7 1 for normal random variables A/"(/i, a 2 ). 

26 



/ dt dxa(t,x) = - I dt dx dx'a(t,x)D(x,x';t,TFR)cr{t,x') (39) 
J t Jt * J t Jt Jt 

and consequently 



dxa(t,x) = - / dxdx'a(t,x)D(x,x';t,TFji)a(t,x') (40) 
2 Jt 



This no arbitrage condition must hold for any Treasury bond maturing at any time x = T. 
Hence, differentiating the above expression with respect to T produces the final version of the drift 
restriction 



/X 
dx 1 'D(x,x';t,TFR)cr(t,x') (41) 

A. 4 Proof of Proposition 

F(to,U,T) = j 'vAe-^ dxf{t "' x) 

= e ~ Sl dx ffo,x) I <pA e - Ito dt J t T t dxa(t,x)+a(t,x)A(t,x) 



= F(t , U, T)e-^o dt & dxa(t ' x) Z[a(t, x)9(x - Q9(T - x)} 
where Z[J] is defined in equation (j37f). Proceeding further leads to 



!F(to t T) = F(to t T)e~^0 dt ^* dX ^ dx ' a ( t > x ) D ( x > x> >*> T FR) a ( t >^') £ k It * dt It* dxdx'cr(t,x)D(x,x';t,T FR )a(t,x') 

= F(t Q ,U,T)exp{Q^t ,U,T)} (42) 
where flp(to, t*,T) is given by equation (|TT|). 
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